Using the uncertainty principle to design simple interactions for targeted self-assembly 
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We present a method that systematically simplifies isotropic interactions designed for targeted 
self-assembly. The uncertainty principle is used to show that an optimal simplification is achieved 
by a combination of heat kernel smoothing and Gaussian screening. We use this method to design 
isotropic interactions for self-assembly of complex lattices and of materials with functional properties. 
The interactions we derive are significantly simpler than those previously published, and it is realistic 
to discuss explicit experimental implementation of the designed self- assembling components. 
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Fabrication of functional materials through directed 
self-assembly is an attractive possibility for many rea- 
sons, such as scalability, robustness, flexibility, and cost 
efficiency. Two main challenges must be addressed before 
this idea becomes useful in practice. First, we must be 
able to determine the necessary characteristics of the self- 
assembling building blocks. There are several approaches 
to this: extensive computational explorations [lffi], prob- 
lem specific heuristics [MS], and, as in this paper, theo- 
retical methods [9HT2]. Second, we must be able to manu- 
facture the components according to the derived design. 
This second point is problematic. So far the theoreti- 
cally designed components rely on interactions that are 
too complicated to be implemented experimentally. 

The gap between theory and experiment is continu- 
ously reduced by the rapid advances on the experimen- 
tal side where it is now possible to fabricate nano- and 
mesoscale particles with complex and tunable interac- 
tions [T3UT6] . While these results are impressive in them- 
selves, the gap between theory and experiments remains 
wide and is unlikely to be closed from the experimental 
side alone. In this paper we report progress from the the- 
ory side, bringing analytically designed systems closer to 
experimental realization. 

In this Letter we introduce a general principle for how 
to systematically simplify isotropic interaction potentials 
that lead to self-assembly of metamaterials with desired 
properties. The method is based on the uncertainty prin- 
ciple, which limits how short-ranged and smooth a poten- 
tial can be without losing the features in reciprocal space 
that guarantee the self-assembly of the target structure. 
A consequence of the uncertainty principle is that heat 
kernel smoothing of the reciprocal potential is optimal in 
the sense that the interaction range decreases at a mini- 
mal loss of detail in the reciprocal (design) space. 

We first illustrate the simplification method by design- 
ing a self-organized disordered material with suppressed 
diffraction around a chosen frequency, i.e., with near in- 
visibility in the chosen frequency range. The second ex- 
ample, our main result, are interaction potentials that 
self-assemble into 2D Kagome and 3D diamond lattices 



and are much simpler than previously published interac- 
tions. The simplifications are significant enough to al- 
low for speculation on possible experimental implemen- 
tations. 

As a model we consider a system of spherically sym- 
metric (colloidal) particles that self-assemble by minimiz- 
ing the potential energy. A configuration is described by 
its density, expressed as a sum of Dirac delta functions 
centered at the particle positions r^, p(r) = J2i <K r ~~ r *)- 
The energy is defined by a pairwise isotropic potential 
V(r) as 

E= J dYdY'V{\Y-Y'\)p{Y)p{Y') = J dkV(k)\p(k)\ 2 

(1) 

where \p(k)\ dk = f^ =k dk\p(k)\ and V(k) is the 
Fourier transform of the (spherically symmetric) poten- 
tial [17] . In [12] it was demonstrated that ([!]) can be used 
to design potentials for directed self-assembly of crystal 
structures by penalizing competing configurations. How- 
ever, direct implementations of such schemes typically 
result in complicated long ranged interactions. 

Simplifying designed interactions. This paper fo- 
cuses on simplification of designed interaction potentials. 
Arguing that a potential is simple requires a definition 
of complexity. In practice such a definition depends on 
what types of interactions can be realized experimentally, 
and is therefore problem specific and not immediately 
useful in a general setting. Still, there are generic types 
of complications that cause problems for experimental 
implementations, such as potentials with many extremal 
points or with abrupt changes in their derivatives. One 
suggested approach for moving towards realizability is to 
consider only convex interaction potentials [9] [18] . Con- 
vexity avoids the problem of multiple maxima and min- 
ima in the potential but instead moves the complexity 
into its derivatives, which is not necessarily easier to 
implement in experiments. 

Rather than trying to find general requirements for re- 
alizability we focus on simplification of a given potential, 
which we argue is a more tractable problem since it is not 
based on an absolute measure of interaction complexity 
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Recall that can be used to design the potential in 
reciprocal space. Combining the reciprocal design princi- 
ple with the potential's physical realization in real space 
shows that the uncertainty principle is applicable. This 
places limits on the simultaneous localization of V(r) and 
V(k). To see in detail how the uncertainty relation works 
in our context, we start by defining a screening of a po- 
tential as a transformation V(r) — >• £(r)V(r), where £ 
is a positive function monotonically decaying away from 
the origin. For technical reasons we also assume that 
£ can be normalized, J drr^{r) 2 = 1. The effect in re- 
ciprocal space of the multiplicative screening is defined 
by a convolution: J*[£(r)F(r)] = Jdk'fdk - k'|)V"(|k|). 
The standard uncertainty relation gives a central result 
that limits the dispersion around zero of £ and its Fourier 
transform £ in d dimensions [T7] . 



var[£ 2 ] var[| 2 ] > 
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The lower bound is achieved exactly when £(r) is a 
Gaussian function, which defines a heat kernel smooth- 
ing in reciprocal space. The dispersion var[£ 2 ] measures 
the degree of screening achieved by the transformation 
V(r) — » £(r)V(r), whereas var[£ 2 ] defines the loss of dis- 
tinction by the corresponding kernel convolution. The 
interpretation of the uncertainty relation is therefore that 
a heat kernel smoothing in the reciprocal space achieves 
minimal loss of distinction for a given degree of screening 
of the potential in real space, and vice versa. 

Based on this argument, we propose the following op- 
eration to simplify a designed potential: 



Cm\ smooth 
V(k) > 



V T (k) = J dk / e- |k - k/|2/r F(|k / |). (3) 



This results in a screening of the potential in real space 



V(r) 



screen 



V T (r) = T~ x \v T ] = e- rr2 /V(r), (4) 



where we ignore the irrelevant normalizing scale factors. 
The role of V(r) and V(k) can also be reversed, i.e., a 
screening in reciprocal space results in a smoothing in 
real space. Both these operations reduce the complexity 
of the interaction potential and they will be used in com- 
bination later in the demonstration with self-assembling 
lattices. The conclusion is that Eqs. ([3|, Q, and their re- 
ciprocal duals define a general method for simplification 
of a designed isotropic interaction. 

We use Gaussian screening and smoothing because it 
is optimal in the sense defined by the lower bound in the 
uncertainty principle. In practice other types of screening 
may be easier to implement, for example an exponential 
damping £(r) = e~ rr like that of the Coloumb interac- 
tion between charged colloids due to free counterions. It 



is straight forward to generalize Eqs. (|3| and Q by con- 
sidering other smoothing kernels and their radial Fourier 
transforms as screening functions. 

Metamaterial with suppressed diffraction. The 
most natural functional property for a self-assembling 
material is designed diffraction properties, since it is de- 
termined immediately by the assembled structure. The 
Fourier transform of a particle configuration (the struc- 
ture factor), tells to what extent diffraction of incident 
light occurs in different directions. In this example 
we demonstrate metamaterials with suppressed structure 
factors within a desired frequency range, self- assembled 
from interactions of varying complexity. Since no scat- 
tering occurs for certain wavelengths in this material, it 
would essentially be invisible at those wavelengths [20 . 

From the expression of the energy in it follows di- 
rectly how to construct a potential that energetically dis- 
favours configurations without suppressed structure fac- 
tor: by letting V{k) be large for wave numbers k in the 
bandwidth we want to penalize. Such potentials have 
previously been used in optimization algorithms to con- 
struct materials with similar properties [21]. The func- 
tional form of the resulting potentials are linear combina- 
tions of a few Bessel functions that oscillate and have slow 
asymptotic decays, V(r) ~ cos(r)/ v /r. They are there- 
fore too complicated and long ranged to be useful in a 
self-assembly context, see Fig. [TJi. This provides a natu- 
ral setting for our method. We can apply the heat kernel 
smoothing in Q to V(k) to construct interactions with 
significantly simpler functional forms. Here the trade-off 
is particularly transparent as the reduced localization in 
V(k) makes the suppressed region of the structure factor 
less distinct. 

To be more specific, consider a particle system in two 
dimensions for which we want to design a potential that 
results in a structure factor suppressed suppressed a spe- 
cific wave length n. Starting with Vo jK (k) consisting of a 
Dirac delta function located at k and applying the heat 
kernel smoothing with parameter r results a potential 
with reciprocal form 



V TtK (k) = e-^+^/o (2k/r) , 



(5) 



where Io is a modified Bessel function and k is the ap- 
proximate location of the maximum (assuming n/^/r is 
large), see Fig. [TJl (inset). From the uncertainty prin- 
ciple it follows that this construction leads to minimal 
range of the potential in real space for a given degree of 
localization around the maximum in reciprocal space. To 
prevent the particles from collapsing into a single aggre- 
gate we follow the idea in [21 and add a low frequency 
penalty to the energy spectrum [22] . 

In Fig. [I] different levels of smoothing are demonstrated 
together with the resulting configurations and their cor- 
responding structure factors. Note that even with very 
strong screening, where the potential shows only a single 
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FIG. 1. Heat kernel smoothing simplifying interactions for suppressed diffraction around a wavelength. A potential is 
constructed using ^ to penalize a certain wavelength k and the inverse Fourier transform gives the corresponding interactions, 
which have long-ranged oscillations (a). By applying the heat kernel smoothing in (|3j, we can successively simplify the 
interactions (b-c) until only a single minimum and maximum remains, (d-f) Low temperature configurations from Monte 
Carlo simulations verifying that the interactions cause self-organization into states with diffraction patterns with pronounced 
frequency gaps (insets show their diffraction patterns). 



over-damped oscillation, the suppression remains signifi- 
cant, as seen in the structure factor. 

Complex lattices. One of the draws of directed 
self-assembly is the ability to design and implement mi- 
crostructure displaying functional behavior, with self- 
assembly of for example chiral catalysts [24] and optical 
crystals [25] achieved to date. Thus much of the theoret- 
ical work in this area has concentrated on solving the in- 
verse statistical-mechanical problem for low-coordinated 
crystals [TT], i.e., designing interactions that lead to self- 
assembly of a target lattice, which is also our focus here. 

Suppose we want to construct an isotropic potential 
that leads to self-assembly of a target lattice. The lat- 
tice can be described in reciprocal space by the structure 
factor |p(k)| 2 , which is positive and has support only at 
discrete points Iq. A basic construction of a potential 
that has the target lattice as a lowest energy state, with 
energy zero, is to let V(k) be a positive function with 
zeros at all points |k^|. This scheme works in principle 
but implemented naively it leads to very complicated and 
long range interactions [12] , illustrated by the gray curves 
in Fig. [2] Starting from the basic design however, it is 
possible to use the heat kernel smoothing and screening 
to produce a much simplified potential that still causes 
assembly of the target lattice. Explicitly we use a con- 



struction on the form 



V T *(r) 



-ak 2 



where V e ^.\ is defined in ([5| [26 and the parameters r 
and <j controls the degree of screening in real and recip- 
rocal space respectively, and e is small. The parameters 
di and |ki| are set so that V ra (k) fulfills the basic design 
principle in [12 , i.e., V T(J (\k\) > and VU(|ki|) = 0, 
with a small perturbation that breaks the energy degen- 
eracy to favor the target lattice. The parameters r and a 
should be chosen as large as possible without losing the 
target lattice as the lowest energy state. 

In Fig. [2] we demonstrate the method on a Kagome 
lattice and a diamond lattice. The main results are two 
maximally simplified potentials that are shown to self- 
assemble into the target lattices in Monte Carlo simu- 
lations. We compare the results to our previously pub- 
lished potential for the Kagome lattice, which was the 
simplest interactions we could achieve with ad hoc tun- 
ing in [T2] . and to interactions that gives a diamond lat- 
tice, derived by Rechtsmann et al. [23 using an iterative 
design scheme (blue curves in Fig. [2| . This indicates how 
effective the systematic approach we suggest is. 
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FIG. 2. The main results: simplified interactions designed for lattice self-assembly. The curves show the interactions from ([6| 
for targeted self-assembly in reciprocal space (inset) and in real space of a Kagome lattice (a) and a diamond lattice (b). The 
gray curves are the basic construction with almost no screening or smoothing while the red curves show the optimized potentials 
after applying the method defined in As a reference we also show previously published potentials (blue curves) causing 
self-assembly of those lattices [T2j [23] . 



Summary and discussion. We have demonstrated 
how the uncertainty principle can be used to design 
relatively simple isotropic potentials for targeted self- 
assembly. Our result also indicate fundamental reasons 
for why an isotropic potential must be relatively compli- 
cated to achieve self-assembly of a complex lattice. 

From the experimental side there is a continuous pro- 
gression of techniques for producing nano- and mesoscale 
components with custom designed interactions and 
shapes. Colloids in particular have been identified as 
suited for tailoring interactions of the type we discuss in 
this article. Effective pair interactions between colloidal 
particles arising from e.g. charge, surface features, elec- 
trolytes or polymers in the solution, or dipole moments 
can be combined and tuned [27j [28]. Complex poten- 
tials with many minima will likely be out of reach for 
the foreseeable future but the simplified potentials we 
present here are much more experimentally accessible. 
For example, the diamond potential of Fig. |2|b) could 
be implemented if a Coloumb repulsion could be com- 
bined with depletion attractions [29 . We believe that 
an implementation of this kind would be an important 
step towards practical realizations of designed functional 
materials. 
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